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9 Abstract 

10 This study analyzes changes in the number of fires detected on forest, grass, and transition 

1 1 lands during the 2002-2009 biomass burning seasons using fire detection data and co-located 

12 land cover classifications from the Moderate Resolution Imaging Spectroradiometer 

13 (MODIS). We find that the total number of detected fires correlates well with MODIS mean 

14 aerosol optical depth (AOD) from year to year, in accord with other studies. However, we 

1 5 also show that the ratio of forest to savanna fires varies substantially from year to year. Forest 

16 fires have trended downward, on average, since the beginning of 2006 despite a modest 

17 increase in 2007. Our study suggests that high particulate matter loading detected in 2007 

1 8 was likely due to a large number of savanna/agricultural fires that year. Finally, we illustrate 

19 that the correlation between annual Brazilian deforestation estimates and MODIS fires is 

20 considerably higher when fires are stratified by MODIS-derived land cover classifications. 

21 

22 1. Introduction 

23 Biomass burning is a large source of particulate matter over South America during the dry 

24 season [Artaxo et al., 1998; Andreae et al, 2002], In addition to causing health problems, 
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particles also perturb the radiative balance, affect plant canopies, and modify cloud properties 

2 [e.g., Andreae et al., 2004]. In the tropics, fires are set to clear forests so that agricultural or 

3 pastoral lands may be developed, but they are also set on existing agricultural or pastoral 

4 lands for nutrient mobilization, pest control, and removal of brush and litter accumulation 

5 [ Crutzen and Andreae, 1990]. A number of studies to date have analyzed fire counts and fire 

6 characteristics over the Amazon [Setzer and Pereira, 1991; Kaufman et al., 1994; Koren et 

7 al., 2007, 2009; Torres et al., 2010], and some have analyzed fire counts over different land 

8 cover types [Eva and Lambin, 2000; Eva and Fritz, 2003; Schroeder et al., 2005; Morton et 

9 al., 2008]. Our study differs from previous studies in that moderate-resolution land cover data 

10 from across the entire Amazon basin are utilized to stratify fires by land cover, for all years 

1 1 between 2002 and 2009. Our study also compares satellite measurements of fires with aerosol 

12 loading, and attributes the source of the aerosol loading to forest and savanna fires 

13 interannual ly. 

14 Studies have found higher correlations between fires and deforestation if land cover 

15 information is employed [Eva and Fritz, 2003]. Forest fires may be used as an indicator of 

16 deforestation in the Amazon [Setzer and Pereira, 1991; Eva and Fritz, 2003], although these 

17 fires are directly associated with some causes of deforestation (i.e. clear-cutting) and not 

18 others (i.e. selective logging). Selective logging, which can increase a forest’s flammability, 

19 may lead to fires that are associated not with deforestation, but rather with forest degradation 

20 [Cochrane and Laurance, 2002; Fearnside, 2005; Asner et al., 2006]. However, since the 

21 majority of Amazonian deforestation is attributed to land clearing for cattle ranching or 

22 agriculture and forest fires are routinely associated with these types of deforestation, a 

23 relationship between forest fires and deforestation is likely to be observed [Setzer and 

24 Pereira, 1991; Morton et al., 2006; Morton et al., 2008]. 
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2 2. Data and Methods 

3 This study explores trends in detected MODIS fires between 2002 and 2009 over the Amazon, 

4 during the biomass burning months of June through November. The region of interest 

5 extends from coast to coast, between 0° and 20° S latitude and between 82° W and 34° W 

6 longitude. The relatively high temporal and spatial resolution of the MODIS sensor, along 

7 with predominantly clear or partially cloudy skies due to persistent high pressure during the 

8 South American dry season, allows for daily detection of fires and aerosols [Giglio et al., 

9 2003 ; Levy et al. , 2007] . 

10 Daily 1-km Collection 5 Level 3 MODIS fire detections (MOD14A1) from the Terra 

1 1 satellite are employed in the analysis [Giglio et al., 2003; Giglio et al., 2006], accessed from 

12 The Land Processes Distributed Active Archive Center (http://lpdaac.usgs.gov). The morning 

13 Terra overpass is used instead of the afternoon Aqua overpass because convective clouds are 

14 less developed in the morning. Fires with a detection confidence lower than “high” (80%) are 

15 excluded from the analysis due to possible false detections in small forest clearings [Morton et 

16 al., 2008; Schroeder et al., 2008]. Fire detections are stratified into three land cover 

17 categories: (1) forest, (2) savanna, and (3) transition between forest and savanna, according to 

18 MODIS-derived land cover categories from the previous year. The Land Processes 

19 Distributed Active Archive Center provides a combined Terra/Aqua yearly land cover product 

20 (MCD12Q1). This product uses decision tree classification algorithms and training data to 

21 assign land cover categories at 500-m resolution [Friedl et al., 2010]. A majority filter is 

22 applied to the 500-m land cover data to upscale the data to the 1-km fire product resolution. 

23 The savanna category contains shrubland, savanna, grassland, cropland, and sparsely 
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1 vegetated and barren land. The transition land cover category is defined as 1-km 2 pixels that 

2 do not have a majority land cover (i.e. 50% split between land cover categories). 

3 Trends in aerosol loading are computed using the MODIS Collection 5 Level 3 aerosol 

4 optical depth (AOD) at 550 nm [Levy et al., 2007]. The monthly mean AOD product from the 

5 Terra satellite (MOD08 M3) is employed to calculate a spatially- and temporally-averaged 

6 AOD over our study region during the biomass burning season. Aerosol data is accessed from 

7 the Level 1 and Atmosphere Archive and Distribution System 

8 (http://ladsweb.nascom.nasa.gov/). 

9 To compare yearly Brazilian fires detected by MODIS to measured Brazilian 

10 deforestation rates, deforestation estimates from the Brazilian Instituto Nacional de Pesquisas 

1 1 Espaciais (INPE) are employed (unpublished data, 2006; 

12 http://www.inpe.br/noticias/arquivos/pdf/prodes2009tabelal.pdf). These estimates are 

13 computed using the high-resolution Landsat sensor, as well as other satellites such as the 

14 China-Brazil Earth Resources Satellite (CBERS), particularly when cloud cover is present. 

15 

1 6 3. Results 

17 Figure 1 shows the interannual trend in daily MODIS fires, summed over all days between 

18 June and November, in our study region between 2002 and 2009. Also plotted is the 

19 interannual trend in average monthly AOD. Previous studies have shown a similar correlation 

20 using AVHRR fires and MODIS AOD [Koren et al., 2007], as well as MODIS fires and OMI 

21 AOD [Torres et al., 2010]. A Student’s two-tailed t-test indicates that total fires and seasonal- 

22 regional mean AOD are correlated at the 99% confidence level, with a correlation coefficient 

23 of 96%. The high correlation suggests that the majority of the aerosol loading is due to 
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1 biomass burning fires, and not other sources such as dust or transported aerosol [Koren et al., 

2 2007]. A map of fire locations for each year is depicted in Figure SI in the auxiliary material. 

3 The sharp decrease in both fires and aerosol loading in 2006 has, in some studies, been 

4 attributed to a tri-national alert service for wildfires enacted in 2006 in the wake of the severe 

5 2005 drought [ Koren et al., 2007, 2009]. Other studies attribute this decline largely to 

6 positive precipitation anomalies that year [Schroeder et al, 2009; Torres et al., 2010]. Figure 

7 1 indicates that a high burning year followed in 2007, suggesting that the tri-national alert 

8 service may not have been the cause of the 2006 decline. Figure 1 shows the basin-wide 

9 interannual variability of fires and smoke, but it cannot explain the increase or decrease in 

10 forest or savanna burning in particular. For that, land cover stratification must be applied. 

11 Figure 2a illustrates yearly fire detections stratified by land cover category. Also 

12 plotted is the ratio of forest fires to savanna fires each year. Within our study region, more 

1 3 total savanna fires are detected than are forest fires or transition fires for each year. Thus, it is 

14 probable that a significant portion of the aerosol loading in our study region is due to savanna 

15 fires, in addition to forest fires. Forest fires generally exhibit a higher particulate matter 

16 emission factor than savanna fires [Janhdll et al., 2010], and also contain more biomass per 

17 area than savanna or other land use types that fall within our “savanna” category, yet 

18 combustion efficiencies are often lower [Ward et al., 1992], If forests were contributing the 

19 bulk of particulate matter each year, we would expect the correlation between forest fires 

20 alone and AOD to be equal to or higher than the correlation between total fires and AOD, 

21 which it is not. The correlation between forest fires and AOD is 0.87 and the correlation 

22 between savanna fires and AOD is 0.85. Both correlations are lower than the correlation 

23 between total fires and AOD, 0.96, suggesting that savanna fires contribute appreciably to 

24 aerosol loading in the region. 
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1 Figure 2a depicts a general decreasing trend in forest fires between 2006 and 2009. 

2 Forest fire detections were only 30% higher in 2007 than in 2006, yet savanna fires were 

3 300% higher in 2007 than in 2006. As a result, the 2007 peak in aerosol loading observed in 

4 Figure 1 is likely due to an anomalously high number of savanna fire detections that year. 

5 Large negative precipitation anomalies occurring over the southeastern Amazon and Cerrado 

6 regions in 2007 match the location of anomalously high numbers of savanna fires that year 

7 (auxiliary Figure SI, Figure S2) [. Janowiak and Xie, 1999; Morton et al., 2008; Torres et al., 

8 2010]. Similarly large negative precipitation anomalies occurred in 2005 (Figure S2), along 

9 with a similarly high peak in AOD (Figure 1), yet the ratio of forest to savanna fires was over 

10 twice as large in 2005 than in 2007, indicating an inherent difference in the type of burning 

1 1 which occurred during these two years. Figure S2 indicates that the drought was more severe 

12 and more extensive in 2005 than 2007, which may partly account for the difference in the 

1 3 forest fire to savanna fire ratio between the two years. 

14 Over the eight year data set, the highest forest fire to savanna fire ratio occurred in 

1 5 2004 (0.62) and the lowest ratio occurred in 2008 (0.24). In 2008, forest and savanna fires 

16 declined by nearly 60% compared to 2007, and declined another 50% between 2008 and 

17 2009. This large decrease in both forest and savanna fire detections resulted in the lower 

18 aerosol loadings observed in 2008 and 2009 (Figure 1). This decline may have resulted from 

19 a combination of moderately positive precipitation anomalies, particularly in 2009, as well as 

20 stricter fire regulations and/or a declining economy; however, it is difficult to conclusively 

2 1 determine the cause of the reduction [ Torres et al. ,2010]. 

22 Figure 2b shows fires stratified by land cover category, normalized by the area of each 

23 land cover category. On average between 2002 and 2009, there were 5,226,760 km of forest, 

24 5,154,306 km 2 of savanna, and 305,025 km 2 of transition land in our study region. For each 
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1 year, the number of fires over transition land per square kilometer exceeded that over both 

2 forest and savanna. This result agrees with previous studies of deforestation in the Amazon - 

3 a large fraction of the land clearing and consequential burning occurs on the edge of the forest 

4 frontier [Skole and Tucker, 1993; Cochrane, 2003]. Increased flammability along forest 

5 edges may also play a role [Cochrane and Laurance, 2002], Between 2002 and 2009, one fire 

6 was detected for every 173 km 2 of forest, 71 km 2 of savanna, and 47 km 2 of transition land. 

7 The methodology of stratifying fires by land cover type also results in higher 

8 correlations between fire detections and deforestation. Figure 3 shows scatter plots and linear 

9 least-squares regressions between MODIS fires summed over each biomass burning season 

10 and INPE estimated deforestation for the country of Brazil. The correlation between total 

1 1 fires and INPE deforestation is relatively low, 0.30 (Figure 3a). However, when only forest 

12 fires are retained, the correlation is much higher, 0.65, which is significant at the 90% 

1 3 confidence level (Figure 3b). 

14 The middle and bottom row of panels in Figure 3 represent the number of fires 

15 occurring on two-or-more and three-or-more days within each 1-km 2 pixel, respectively, for 

1 6 each biomass burning season. These higher frequency fires are likely to be associated with 

17 deforestation, as deforested wood is often piled and rebumed on several occasions following 

18 the initial bum [Morton et al., 2008]. The stratified correlations in Figures 3e and 3f, 0.77 and 

19 0.76, are significant at the 97% confidence level. Yet, the difference between the stratified 

20 and non-stratified correlations in this second row is less than between Figure 3a and Figure 

21 3b/c. 

22 The bottom row of panels in Figure 3, which represents three-or-more fire detections 

23 per pixel, shows nearly no difference in the correlation coefficients between the stratified and 

24 non-stratified cases, yet correlations in all panels are significant at the 99% level. 
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1 Stratification improves the correlation to a lesser extent in the case of persistent fires because 

2 persistent fires are often classified as forest or transition fires. Summed fires between 2002 

3 and 2009 indicate that 53% (65%) of two-detection (three-detection) persistent fires occur on 

4 forest or transition lands, compared with 33% of all fires. This method of correlating fires to 

5 deforestation performs well in areas with high rates of clear-cutting deforestation, such as the 

6 states of Mato Grosso or Para [Fearnside, 2005; Morton et al., 2006]. One should take 

7 caution extending this methodology to other regions where deforestation and biomass burning 

8 are not as closely linked, and where humid conditions prevent the widespread use of fires for 

9 land management [Eva and Fritz, 2003]. 

10 

1 1 4. Conclusions 

12 This work demonstrates a statistically significant correlation between yearly MODIS AOD 

13 and MODIS fire detections during the Amazon biomass burning seasons between 2002 and 

14 2009. Stratifying fire detections by MODIS land cover type provides further information into 

1 5 the source of the aerosol loading. A decreasing trend in forest fires is observed between 2006 

1 6 and 2009, despite a modest increase in 2007. Results also indicate that high regional aerosol 

17 loading detected in 2007 was likely caused by a high number of savanna/agricultural fires that 

1 8 year, in part due to negative precipitation anomalies in the southeastern Amazon and Cerrado 

19 regions. We also find that the number of fires occurring in transition regions, regions on the 

20 edge of forest and savanna lands, is higher than both forest and savanna fires on a per area 

21 land cover basis. Finally, we illustrate that the correlation between detected fires and annual 

22 Brazilian deforestation estimates from high-resolution sensors becomes greater if fires are first 

23 stratified by land cover type; however, when only persistent fires are employed, the 
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1 correlation does not improve with stratification since the majority of persistent fires occur 

2 over forest and transition lands. 

3 We acknowledge the various uncertainties associated with this analysis. A number of 

4 fires may not be detected due to atmospheric conditions, satellite overpass time, or fire size 

5 and energy. False fire detections may be included in the data, however, only high-confidence 

6 fires were selected to minimize this error. Misclassified land cover data, especially near 

7 regions of forest conversion, may also introduce potential uncertainty. Yet, this research 

8 indicates that segregating fires by moderate-resolution land cover data provides a simple 

9 method of gaining additional information about the source of regional aerosol loading and 

1 0 about regional deforestation on interannual timescales. 
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2 Figure 1 : Interannual trend of seasonally-averaged MODIS Terra AOD and interannual trend 

3 in total MODIS Terra daily fire detections over the Amazon biomass burning season (Jun - 

4 Nov), for the years between 2002 and 2009. 
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2 Figure 2: (a) Interannual trend of MODIS fire detections during the Amazon biomass burning 

3 season, stratified by MODIS land cover categories from the previous year, between 2002 and 

4 2009. Also plotted is the ratio of forest fires to savanna fires for each year, (b) Interannual 

5 trend of MODIS fire detections stratified by land cover category, per total area of each land 

6 cover category in the study domain. 
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Figure 3: (a) Scatter plot and linear least-squares regression between MODIS detected fires 
over Brazil and Instituto Nacional de Pesquisas Espaciais (INPE) Brazilian deforestation 
estimates for the years between 2002 and 2009. (b) Same as (a) but only including fires over 
the forest land cover category, (c) Same as (a) but only including fires over the forest and 
transition land cover categories, (d)-(f) Scatter plots between two-day persistent fires, defined 
as fires with two or more detections occurring over the same 1-km 2 pixel in a single biomass 
burning season, and INPE deforestation, (g)-(i) Scatter plots between three-day persistent 
fires, defined as fires with three or more detections occurring over the same 1-km pixel, and 
INPE deforestation. The summation of fires between 2002 and 2009 (N), the Pearson’s 
correlation (R), and the non-directional p-value (P) are also included in each panel. 
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